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O ! Abstract 

m 

We address the question of whether or not fermions with twisted periodicity 
condition suppress the semiclassical decay of x Kaluza-Klein vacuum. We 
consider a toy (l+l)-dimensional model with twisted fermions in cigar-shaped Eu- 
clidean background geometry and calculate the fermion determinant. We find that 
(—j ■ contrary to expectations, the determinant is finite. We consider this as an indication 

that twisted fermions do not stabilize the Kaluza-Klein vacuum. 

^ , Keywords: Kaluza-Klein theory, semiclassical vacuum decay, twisted fermions, 2- 

' dimensional model. 

oo ■ 
in 

\f\ • 1 Introduction 
o ■ 

^ ' It is known that x Kaluza-Klein vacuum is unstable fl]. The vacuum decay 

proceeds through the Euclidean bounce, whose metric is 



><\ ds' = -^ + edn^ + R'(i-^)de\ (1) 



where R < ^ < oo, < 9 < 2t^ and dVt^ is the metric of unit 3-sphere. At large ^, the 
geometry of this solution is i?^ x S"^, whereas as ^ — )■ i?, the geometry approaches R"^ x S^. 

2 

The latter property is seen by performing the change of variables ^ = R + which gives 
near r = 

ds^ = dr^ + r^de^ + R^d^'' . (2) 

Upon continuation to the space-time of Minkowskian signature, the bounce ([T]) describes 
the decay of x into nothing [T] . The decay of the Kaluza-Klein vacuum and similar 
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processes are of interest from several viewpoints, and various versions of this phenomenon 
have been extensively discussed in literature pHT5]. 

It has been argued in Ref. pTj that x Kaluza-Klein vacuum is stabilized in a 
theory containing fermions with twisted periodicity condition in the compact coordinate, 

^(e) = e-2^-'?^(^ + 27t) . (3) 

The argument is that in the background of the bounce ([T]), this condition makes perfect 
sense away from ^ = i?, but for generic q it appears singular at ^ = i?. The way to 
see the effect of fermions on the Kaluza-Klein vacuum decay would be to calculate the 
ratio of fermion determinants in the background ([T]) and in the background of the original 
Kaluza-Klein metric, both with periodicity condition (|3j): 

det 

A(g) = — — , (4) 

where V is the Dirac operator and index q refers to the condition ([3]). If this ratio 
vanishes, the vacuum is stabilized indeed. 

In this paper we consider a toy model for this situation. Namely, instead of studying 
the five-dimensional theory with the metric ([1]), we discuss two-dimensional theory with 
the metric describing a sigar-shaped space: 



?2 



ds' = + R'\^1- ^ J de' (5) 

In other words, we disregard the part of the geometry. Likewise, instead of the 
five- dimensional Kaluza-Klein vacuum we consider (1 -|- l)-dimensional theory with the 
spatial dimension compactified to a circle, whose Euclidean counterpart is described by 
the metric 

rfsL = d^"^ + R^de'^ . (6) 

The metrics ^ and ([HD coincide as ^ — i- 00, but the geometry ^ has a smooth end at 
C, = R, just like in the case of the five- dimensional bounce. We impose the periodicity 
condition ([3]) on the two-dimensional fermions, which again makes sense away from ^ = 
but appears singular at ^ = i?. Thus, we argue that our toy model captures the main 
features of the five-dimensional Kaluza-Klein theory, which are relevant for the vacuum 
decay. Our purpose is to calculate the ratio (jl]) in this toy model. In our calculations we 
consider the interval 

i? < e < T , (7) 

where T is sent to infinity in the very end of the calculation. We do this for the deter- 
minants in the backgrounds of both "bounce" metric (jSj) and "vacuum" one ([6]). The 
choice of one and the same IR cutoff T in these two cases mimics the five-dimensional 
situation, where ^ is unambiguously defined as the radius of the 3-sphere both for the 
bounce solution and Kaluza-Klein metric. 

To simplify things further, we consider massless fermions. The motivation is that a 
pathology, if any, in the fermion behavior in the background metric ([5]) would emerge 
from a small vicinity of the point ^ = R, while the short-distance properties of fermions 
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should not be sensitive to their mass. The advantage is that we can utihze conformal 
invariance of massless two-dimensional fermions (modulo conformal anomaly, which is 
local, and therefore independent of q). Indeed, our metric ([5]), as any other 2D metric, 
is conformally fiat, 

ds^ = ^''{r)[dr^ + r^de^], (8) 

where 



e « 

$ = R- — — e « 
2e 



(9) 
(10) 
we 



Thus, instead of calculating the fermionic determinant in the background metric ([5]) 
are going to perform the calculation on a plane, still with the periodicity condition 

It is worth noting that the calculation of det 2)^"""'=^ is equivalent to the calculation 
of the determinant for conventional (periodic) massless fermions of charge {—q) in the 
background of an instanton in the two-dimensional Abelian Higgs model, in the limit of 
vanishing instanton size. If there is no interaction of fermions with the scalar field, then 
the fermionic part of the Lagrangian in the latter model has the form 



IV. 



iip-f^{df, - iqeA^)i) , 

while in the limit of vanishing instanton size, the field of the instanton — Abrikosov- 
Nielsen-Olesen vortex [T6l[l7] — has the Aharonov-Bohm form, = €'^8^6. The field 
ip obeys the periodicity condition ip{6 + 2n) = il){9), so the change of variables ip — )• e^'^^ip 
reduces the problem to the calculation of the fermion determinant on R^ without the 
gauge field background, but with the twisted condition 

Determinants in the instanton background are well studied in (1 + 1) dimensional 
models. For scalar and vector fields the calculation was performed, for example, in Ref. 
^18j . and for fermions in Refs. [19l[20]. In Refs. [2T1[22], the calculations were performed 
for chiral fermions of half-integer charge, coupled to the scalar field. Interestingly, the 
determinants of fermions with half- integer charge do not show any pathology [2T1I22]- 

Somewhat surprisingly, in this paper we show that similar result holds in our problem: 
the fermion determinant in the background of our toy "bounce" is finite for arbitrarjl^ 
q 7^ ±1/2. We consider this as an indication that twisted fermions do not, in fact, 
stabilize the x S*^ Kaluza-Klein vacuum. 



2 Fermion determinants 

It is convenient to split the quantity of interest, the logarithm of the ratio (jl]), as 
follows: 



InA(g) = In 



"det V\°'''"'^' 




" det 'j^'"'^'^''^''^^ ' 




" det V^g~ ' 




-In 




+ In 


det V\°^^^^ 


det V™"'_ 


det V^gl~ 



(12) 



Antipcriodic fermions {q = ±1/2) are special, see below, and are not studied in this paper. 
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According to the above discussion, we are going to calculate det on a plane, and 

det T''"^'^""™' on a cylinder ([6]). For q = 0, the fermion behavior in the background metric 
(|5]) is manifestly healthy, so the last term must be finite. The explicit demonstration of 
the latter fact is somewhat subtle; we give the corresponding analysis in Appendix. 

Let us concentrate on the first two terms in the right hand side of eq. f[T^ . Let us note 
that A(g) is periodic in q with period 1, so we can consider, without loss of generality, 
the range —1/2 < q < 1/2. Furthermore, due to C-invariance, A(g) is symmetric under 
q — )■ —q. Therefore, it is sufficient to study the theory with < g < 1/2. We are going 
to perform our calculation for 

< g < i . (13) 

The case g = 1/2 is subtle for reasons that will become clear later, and we leave it for 
the future. 

2.1 Vacuum background 

We begin with the vacuum-vacuum term. We recall that 

In [det p^"^""'"] = -EgT , (14) 

where T is the normalization time and Eg is the Casimir energy of fermions with twisted 
periodicity condition in the x theory. This Casimir energy was calculated in Ref . |2] , 
but the resulting expression there is somewhat implicit. So, we redo the calculation here. 
We write the Casimir energy as a sum over energies of the Dirac sea levels, 

oo 

n=— oo 

where n is the angular spectral number and aj„ = \n + q\/R. We regularize this sum by 
multiplying each term by a factor exp{—e\uJn{q)\) , where e is a small parameter sent to 
zero in the end of the calculation, and obtain 



Eg- Eo = - — 



(16) 



With our convention f ll3p . this sum is straightforwardly evaluated, and in the limit e — )■ 
we obtain 



In 



det 2^"'*'^""™ 



-{Eg - Eo)T = ~{q' - q) (17) 



We have compared numerically this simple result with that of Ref. [2] and found excellent 
agreement. 

2.2 Bounce background 

Let us now turn to the bounce-bounce term in Eq. 0121) . which involves the ratio of 
determinants on a disc r < a, where a is the IR cutoff in terms of the coordinate r, cf. 
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Eq. ([7]). As the operator D^°^^'^^ is anti-Hermitean, its eigenvalues are purely imaginary: 

In det I)^™ = ^ , (18) 



Lm 



where / and m are the radial and angular spectral numbers, respectively, and A^^^ are the 
eigenvalues of the operator p^o""-'^^ 



•\bounce„ 



2^'™^/,™ = Yd.Am = l\mi^l,m • (19) 

We take Euclidean 7-matrices and spinors in the following form: 

and impose the Dirichlet boundary conditions in the radial coordinate, 

0(a,^) = O. (20) 
We require that the eigenfunctions be square-integrabl^. 



The component x is related to as follows 

1 



X = r^(5o0 + ^5i0), (21) 



while obeys 

The eigenfunctions are 

where 



A0 = -AV (22) 
V(^,^)=e™Vn^(r), (23) 



m = n + q (24) 
and n is integer. 4>i,m{r) obeys the Bessel equation 

r'<Pi,m + r<P'Lm + [ALr2 - m^]^,,^ = . (25) 



^bounce 



Obviously, the eigenvalues come in pairs, ±A;_m, where A^^^ > 0. Note, that V 
has no zero mode: even though for A = both and x can be square-integrable, it is 
impossible to satisfy the boundary condition fl20l) . So, we have 



ln[detP^™] = J2 ^^^Im ■ 

The square-integrable solutions are (pi^m = Jm{rXm) = Jn+q{i^Xn+g) for n > and (pi^m = 
J_m(rA_m) = J-ri_g(rA_„_g) for n < —1. It is convenient to rewrite the latter function 



^'This requirement ensures, as usual, that the expansion of the fermion field in the path integral, 
Tp — ai,m'4^i,m, where a;.,„ is the set of the integration variables, yields the diagonal fermion action 
with finite coefficients, Sp — J d?x i/jt2?'""'"ce^ _ m^^i,™-^] m°^i,m ^ d?x ipj ^ipim- Note that the 
operator 2?^'°""'='= is indeed anti-Hermitian on square-integrable solutions obeying (I^Ul) . i.e., that the 
boundary terms appearing when integrating (J cPx 'ip\T>^°^^'^'^'ip2)* by parts, vanish. 
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as J„'_g(rA„'_q), where n' > 1; in what follows rriq denotes either n + q or n' — q. It is 
straightforward to see that for these solutions, the components x given by Eq. fl2T]) are 
square integrable as well. 

Let us point out a subtlety of the antiperiodic fermions, q = 1/2. In that case the 
solutions to Eq. f lT^ with n = —1 are not square integrable: either (p ot x behaves as 
r~^/^ near the origin. We leave the analysis of this case for the future, and proceed with 
the model with < g < 1/2. 

2.2.1 Sum over radial spectral numbers 

The radial spectrum of eigenvalues is determined by the following equations: 

Jn+q{aM,n+q) = ] 72 > (26) 
Jn'-q{a\l,n'-q) = ] u' > I (27) 

Let us evaluate the sum over radial spectral numbers / for a given n using (^-function 
method, 

where 

oo 

Zm{s) = - Uis) = [XiX - Ai:) (28) 

1=1 

The direct calculation of this sum is complicated by the fact that the eigenvalues are 
not known analytically. To this end, a convenient tool is the Gelfand-Yaglom formalism 
[23ti26] . which was successfully applied for calculations of other functional determinants. 

Let us briefly describe the method. If we know a function F{z) which has zeros at 
desired eigenvalues, which are assumed to be positive, then its logarithmic derivative has 
simple poles with residues equal to 1 at those eigenvaluse, and we can write ^-function 
as follows: 

C(s) = J]Res 
I 

Assuming that F{z) is non-singular anywhere except possibly for real negative semi-axis 
and z = oo, this can be transformed into a contour integral, see Fig. [H 

(is) = C^{s) + Cn{s) , (30) 

where 

Us) = 7^ [ dz z-^ ^ In F{z), (31) 
zTTi J az 

7 

n 

The contour 7 surrounds the negative real semi-axis, and f2 is a large circle. The inte- 
gration runs clockwise in (I3T]) and counter-clockwise in (l32l) . 



'-^\nF{z) 

az 



(29) 
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Im z 




The first trial in our case would be Fi^z) = Jrn{0'Z), which has zeros ai z = Xi^m- 
However, we actually need zeros at z = Xf^. Also, we have to avoid a zero at z = 0. The 
function that satisfies these requirements is 



F(z) 



{a^/z) 



Let us first calculate the integral over the large circle. 



2m 



dz z 



d_ 

dz 



In 



{a^/z) 



L^mo 



d 



dz 



(33) 



(34) 



To estimate the behavior of the term containing Bessel functions, we make use of the 
"approximation by tangents" [27] at large \z\ and possibly large index m. 



^miz) 



Try] 
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+ 



128(^2 
5m^ 



cos I yz^ 



m TT 



— m arccos ■ 



sm I Vz^ 



— m arccos ■ 



z 

m TT 
7 ~ 4 



+ 



(35) 



and find that the first term in the integrand in (IMI) behaves as 1 2; | at complex infinity. 
Therefore, its contribution to the integral vanishes as s — )■ 0. The contribution of 
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the second term is straightforwardly evaluated and gives 



ZmAO) = ^ 



where L is the radius of the large circle. 

Let us turn to the remaining part, Zm^^(s). We recall that 



(36) 



Jm{e^x) = e^'~ lm{x); argx = 



(37) 



and get 



•^ni,7 ("S) 



Sm ITS 



TC 



dx X ( — In 

ax 



ImAay/x) 



^moi^y/^) _ dx 



d 



+ In [a"^0-m,^(mo-m,)/2j 



We see that Z'^ .^{0) is an integral of total derivative. The contribution due to the upper 
limit X = L is proportional to InL and exactly cancels out the contribution fl36|) . Taking 
into account that In [im^^ay/x) / lmo{ay/x)~\ vanishes as x — )■ oo, we finally get 



lim In 



In 



m,-mo r(mo + 1) 

2) T(m, + 1) 



(39) 



2.2.2 Sum over angular spectral numbers 



We are now in a position to calculate the ratio of determinants. We regularize the 
sum over angular spectral numbers by multiplying the determinants by exp (ie | | ) , cf. 
Eq. f|T6|) . and write explicitly 



In [det Vl~] exp{ie\de\)- In [det exp{i6\de\) 



n=l 



a\n+q-l 



n'=l 



2, 



r-\n + q) 



^-e{n+q-l) ^ 



-(n-1) 



Tin) 



-e(n-l) 



+ 



T-\n'-q + l) 



^-e(n'-q) ^ 



2. 



r(n' + 1) 



(40) 



where we have shifted the argument in the first sum, n — )■ (n — 1). The right hand side 
of Eq. fHOl) contains a-dependent and a-independent parts. 



ln[det pj™] exp(2e|ae|)-ln[det p^™] exp{ie\de\) = ^ A(ri, g)+5]] ^(n, g, a) , (41) 
where 



n=l 



n=l 



n=l 



n=l 



In r(n)e-"'^ + In r(n + l)e-"("+^) 

- lnr(n + g)e"("+'?) - lnr(n + 1 - g)e"("+^-'?) 



(42) 
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and 

oo oo 

Hn, g, a) = In I ■ ^ [(n + g - i)e-^(-+9-i) - - 1)6" 

+ {n- g)e-"("-'?) 



-e(n-l) 



n=l 



n=l 



ne 



(43) 



The a-dependent sum, Eq. ( 143|) . coincides with that of Section [2?T| and we immediately 
obtain that in the hmit e — )■ this term is 



(g - g2) In ■ 



(44) 



n=l 



We now recall that r = a is the cutoff radius on a plane, which is related to the cutoff 
^ = T by Eq. i.e., ln(a/2) = T/R. Thus, the contribution (jHj) cancels out the second 
term in Eq. (|T2l) . which is given by Eq. (|T7|) . As could have been anticipated, the result 
for InA(g) is infrared finite. 

The quantity of interest, InA(g), is thus given entirely by the a-independent sum fl42p 
(modulo the last, g- independent term in Eq. f|T2|) ). To extract its part which is potentially 
divergent in the limit e — )■ 0, we make use of the Stirling approximation, 

Inr(n) = [\nT{n)]g^ + 0{n-^) , 

where 

[Inr(n)]^, 

Hereafter the subscript St denotes the quantities calculated within the Stirling approxi- 
mation. Since A(n) — Ast{n) decreases as at large n, potentially dangerous are the 
sums involving the Stirling approximation of Inr(n). These boil down to 



n\nn — n In n H — In 27r H 

2 2 Un 



(45) 



n=l 
oo 

n=l 



''^"'n In n 



E 

n=l 



--^ Li^(e^ 
as 

Tnn = — — Lis(e~ 
ds ^ 

p-en 

= -ln(l-e-^) , 

n 



s=0 



(46) 
(47) 
(48) 



where ljis{z) is polylogarithm. We find that the terms that diverge as e — ?■ actually 
cancel out due to the identity 



lim ( e^— Ligie 

e^O V ds ^ 



ds 



+ ln(l 



s=0 
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(49) 



where 7 is the Euler-Mascheroni constant. Thus, the ratio of fermion determinants, A(g) 
is finite, which is our main result. The explicit expression is 



In 



"det 




■ f^Q^ 'j^vacuum " 




-In 




det 


det 



^ 00 

7(g - g') + - (27 - 1 + ^(1 + g) + ^{2 - g)) + 5^(A(n, g) - Ast{n, g)) , (50) 



n=l 
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Figure 2: Logarithm of the fermion determinant as function of the twist q. 



modulo finite and g-independent constant. The leading contribution here comes from the 
analytical part. This function is shown in Fig. |21 

To summarize, we have found that all would-be divergences cancel out, and the de- 
terminant of twisted fermions is finite in the background of our "bounce" . Were similar 
situation inherent in the five- dimensional theory, the Kaluza-Klein vacuum would be 
unstable even in the presence of twisted fermions. 
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Appendix 

Let us show that the last term in Eq. f|T2|) . corresponding to periodic fermions, is 
finite. We begin with the vacuum term and again use Eq. ( IT^ . now with g = 0. We write 

2 °° 

n=l 

where e is the regularization parameter. It is straightforward to obtain 

_ 2 1 
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The interpretation of the first, divergent term is that it corresponds to the vacuum energy 
density on infinite hne, which should be renormahzed away. Indeed, physically meaningful 
is the cutoff in energy, exp(— a;„/A). Since Un = \n\/R, we identify e = {AR)~^, and the 
first term becomes E = —2h?R. The corresponding energy density E/R is independent 
of i?, so it is indeed the energy density on a line. So, the Casimir energy of periodic 
massless fermions equals 

^ ~ m' 

Let us consider now the bounce term. We introduce the notation a;(r; R) = In $(r; R), 
and consider the variation of the determinant under the change of the value of R. The 
conformal anomaly gives 



^fllndetP^""""" = -— I d^x 6ruj -Du = — I d^x d^Snco ■ d^u 



If 1 
— / d^x 5ruj -Du = — 

OTT J OTT 

1 r 

/ d9 a6R[ijj{r = a)]9rw(r = a) . (51) 

OTT J 

A subtlety here is that the boundary term does not vanish and, in fact, it is important 
for the cancellation of the infrared divergence. At large r we have u = — ln(r/i?), 
drUj{a) = —1/a and 

T 

oj{r = a) = — lna + lni? = — — + finite , 

R 

where we recalled that In a/2 = T/R and omitted the terms which are finite in the limit 
T — 7- oo. We immediately obtain 

In det = / d^x d^u ■ d^u - ^ + finite . 

l/TT J 6R 

This expression shows that there is no UV divergence in In det V^'^'^'''', as expected. The 
IR divergent part of the first term here equals (1/6) In a = T/{6R). Thus, 

In det V^"'"'"^^ = + finite . 

6R 

Its IR divergent part is equal to In det 2)q"™"™ = —ET, so the last term in Eq. ( IT2|) is 
indeed finite. 
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